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ABSTRACT 

We analyze the properties of dark matter hales in the cold-plus-warm dark matter cos- 
mologies (CWDM). We study their dependence on the fraction and velocity dispersion 
of the warm particle, keeping the free-streaming scale fixed. To this end we consider 
three models with the same free-streaming: (1) a mixture of 90% of CDM and 10% of 
WDM with the mass 1 keV; (2) a mixture of 50% of CDM and 50% of WDM with the 
mass 5 keV; and (3) pure WDM with the mass 10 keV. "Warm" particles have rescaled 
Fermi-Dirac spectrum of primordial velocities (as non- resonantly produced sterile neu- 
trinos would have). We compare the properties of halos among these models and with 
ACDM with the same cosmological parameters. We demonstrate, that although these 
models have the same free-streaming length and the suppression of matter spectra are 
similar at scales probed by the Lyman-a forest (comoving wave-numbers fc < 3 — 5 
h/Mpc), the resulting properties of halos with masses below ~ IO^^Mq are different 
due to the different behaviour of matter power spectra at smaller scales. In particular, 
we find that while the number of galaxies remains the same as in ACDM case, their 
density profiles become much less concentrated. Our results imply that a single pa- 
rameter (e.g. free streaming length) description of these models is not enough to fully 
capture their effects on the structure formation process. 
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1 INTRODUCTION 

It is usually said that cosmological data favour Cold 
Dark Matter. The cosmological "concordance model" is 
therefore often called ACDM. However, a more precise 
statement should be that hot dark matter particles (i.e. 
the particles that became non-relativistic only around 
recombination time, such as e.g . the ordinary neutrinos) 
are ruled out (|Davis et al.l Il985l ) . The difference between 
cold and warm DM particles (the former being always 
non-relativistic and the latter becoming non-relativistic 
deeply in the radiation-dominated epoch) would show up 
at approximately galactic scales and it is only recently 
that such small scale effects are sta rting to be resolv- 
able bo th theoretically (se e e.g. iMaccio &: Fontanot 
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ulations require significantly larger number of particles to 
reso lve the same scales as compared with the CDM case (see 
e.g. IWang fc Whit3l2007l : iLovell et al.ll201ll ). Additionally, 
at sub-Mpc scales baryonic physics ca n hide (or mimic) th e 
WDM suppre ssion of power (see e.g. Benson et al.l (|2002l ): 
iBuUock et aD (|200d ): ISemboloni et all (|201ll ). which makes 
the analysis of small-scale data challenging. 

The tiny (from the cosmological point of view) differ- 
ence between cold and warm dark matter is however of cru- 
cial importance for particle physics, as it means a huge dif- 
ference in the properties of corresponding particles and may 
eventually provide a clue on the structure of a fundamental 
theory of particles and interactions. 

Historically the first WDM models were thermal relics 
- particles that were in equilibrium in the early Uni- 
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verse and f roze-o ut, being relativistic (| Colin et al.l I2OO0I : 
iBode et al.l I2OOII ). Such particles had thermal primor- 
dial velocity spectrum and strong cuto ff - like suppression 
of the power spectra (|Bode et all I2001I : IViel et all l2005l : 
iBovarskv et al.|[2009al ) at scales below few Mpc. Such models 
are characterized by only one scale — position of the cut- 
off in the power spectrum, related to their free-streaming 
horizon. One possible tool to probe the growth of struc- 
tures of (sub)Mpc scales is the Lyman-a forest method - 
studies of statistics of absorption features in the spectra of 



distant quasars. The Lyman -a forest data ([Hansen et al 
2002; Vicl ct al. 2005, 200d|; ISeliak et al.l l2006l : IViel eTal 
20081 : IBovarskv et al.l |2009a|) puts such strong constraints 



at their free-streaming length that "thermal relics" WDM 
models, compatible with Lyman-a bounds, produce essen- 
tially no observable changes in the Galac tic structures (c.f . 
Strigari et al. 20_06; CoUn ct al. 2008; Bov arskv et al.ll2009al : 
de Narav et aril2009i : ISchneider et al.ll201ll v 



However, particle physics motivated WDM candidates 
can be produced in the early Universe in non-thermal 
ways, may have significant non-zero primordial velocities 
in the radiation-dominated epoch and non-equilibrium v e- 
locity spectra (for review see e .g. IBovarskv etal] l2009bl l3: 
iTaoso et all l2008l : iFend I2OIOI '). In many models (e.g. 
sterile neutrinos, gravitino, axino) the same DM par- 
ticles can be produced via two co-existing mechanisms 
and therefore generically primordial velocity spectra have 
"colder" and "warmer" components. Such models can be 
called mixe d or "cold plus warm " dark matter models 
(CW DM) (|Bovarskv et al.l l2009al . see also iPalazzo et al.l 
I2OO7I ). Qualitatively, structures form in these models in a 
bottom- up fashion (similar to CDM). The way the scales 
are suppressed in CWDM models is more complicated (and 
in general less severe for the same masses of WDM par- 
ticles), as co mparable with pur e warm DM models. The 
first results of iLovell et al.l (|201ll ) demonstrate that the res- 
onantly produced sterile neutr ino DM models , compa tible 
with the Lyman-a bounds of IBovarskv" et al.l (|2009bl ). do 
change the number of substructure of a Galaxy-size halo 
and their properties. The discrepancy between the number 
of observed substructures with small masses a nd those pre- 
dicted by ACDM model s (first pointed out in iKlvpin et al.l 
I1999I : iMoore et al.lll999l ) can simply mean that these sub- 
structures did not confine ga, s and a r e therefore completel y 
dark (see e g. iBullock et al.l (I2OO0I): iBenson et all j2002l ): 
ISomervilld (|2002l '): iMaccio et al.l (|2010l ')'). This is not true 
for larger objects. In particular, CDM numerical simulations 
invariably predict several satellites "too big" to be masked 
by galax y formation processes, in co ntradiction with obser- 
vations (|Bovlan-Kolchin et al.|[201ll ). Sterile neutrino DM 
of the mini mal neutrino extension of the Standard Model 
( "the i/M SM) lAsaka fc Shaposhnikovl (|2005l ') ; IBovarskv et all 
l|2009cl ). with its non-trivial velocity dispersio n, turns out 
to be "warm enough" to amend these issues (|Lovell et al.l 
I2OIII ) a nd "cold enough" to be in agreement with Lyman-a 
bounds (|Bovarskv et al.ll2009bl ). 

In this paper we study the structures of halos in the 
CWDM models. To this end we pick two CWDM mod- 
els, com patible with the Lyman-a data of IBovarskv et al] 
l|2009al ) and having the same free-streaming of the WDM 
particles. We demonstrate that the properties of halos differ 
in these models (and differ from both pure CDM and pure 
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Figure 1. Ratio of the power spectra of the three models used in 
this work to that of ACDM (for the same values of cosmological 
parameters). The blue short-dashed curve is a mixture of 90% of 
cold dark matter and 10% warm dark matter with rescaled Fermi- 
Dirac spectrum and mass mwDM = 1 keV, (model ml); the red 
solid line is a mixture 50% of cold and warm particles with the 
mass mwDM = 5 keV (model m5) and the green dashed-dotted 
line (model wdm) is a 100% WDM model with the mass of DM 
particle mwDM = 10 ke V. All three models are compatible with 
the Lyman-a analysis of lBovarskv et al] ((2009^. 



WDM model with the same free-streaming), meaning that 
they are not determined by the free-streaming alone. 

The paper is organized as follows. We discuss the choice 
of parameters of our DM models and initial conditions for 
simulations in Section[2] The suite of simulations is discussed 
in Sec. [3] and the main results in Section |4l We discuss our 
results in Section (5] 



2 MODEL SELECTION AND INITIAL 
CONDITIONS 

For our simulations we selected two CWDM and one WDM 
models. Our first CWDM model (model ml in what follows) 
is a mixture of 90% of cold dark matter and 10% of warm 
dark matter with the mass mwdm = 1 keV, while the second 
is a 50% cold and warm mixture with the mass mwdm = 
5 keV (model m5 hereafter) . The WDM model has 
10 keV. In all cases the WDM components have the rescaled 
Fermi-Dirac spectrum of primordial velocities 
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exp 
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where m^dm is the WDM mass, T^{z) is the temperature 
of cosmic neutrino background, evolving with redshift as 
T„{z) = Tvo{l + z), Tvo = 1.9 K, and the constant x is 
determined by the requirement to provide a given fraction, 
Jwdm, of the total dark matter density: 
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(in the units where c = h = ks = !)■ Non-resonantly pro- 
duced sterile ne utrinos would have such a ph ase-space dis- 
tribu tion ^ (jPodelson fc Widro"w|[l994, : Dolgov fc HansenI 
|2002| ). One can relate the para meters of such W DM models 
and those of "thermal relics " (|Bode et al.ll200ll ). using e.g. 
formulas in IVielet al.l (|2005l ). In order to investigate the de- 
pendence of the properties of the cosmologies on parameters 
other than free-streaming (i.e. WDM mass and F^^m com- 
ponent), we have chosen the models in such a way that their 
free- streaming horizon is the same. Both mo dels are compat- 
ible wi th the CWDM Lyman-a analysis of iBovarskv et all 
|2009al ). 

We computed the linear power spectrum at redshift 
Zini ~ 30 of matter density perturbations Pin ijk, Zini) for 
those models. The standard software (i.e. CAMB. [Lewis et al] 
[2000) is not immediately appropriate for this purpose, as it 
only treats massive neutrinos with a Fermi-Dirac primor- 
dial distribution. To adapt it to the problem at hand, we 
modified CAMB so that it could take arbitrary spectra as in- 
put data files. We analyzed the spacing in momentum space 
needed in order to obtain precise enough results, and imple- 
mented explicit computations of distribution momenta in 
CAMB. We cross-checked our results by modifying another 
linea r Boltzmann solver - cmbfast (jSeliak fc Zaldarriagal 
Il996t ). implementing a treatment of massive neutrinos with 
arbitr ary analytic distribution fu nction and with the CLASS 
code (|Lesgourgues fc Tramllioill ) where such an option is 
realized. 

Figure [1] shows the ratio of the power spectrum in our 
CWDM models with respect to the ACDM model with 
the same cosmological parameters. The effect of the free- 
streaming is clearly visible on scales smaller than fc « 1 
h/Mpc. It is important to notice that although the two mod- 
els formally have the same free-streaming length, since the 
product of the WDM mass and its abundance is constant, 
the suppression of power on small scales is different in the 
two models and larger for the m5 model, due to its more 
abbundant warm dark matter component. 

Based on the linear power spectrum Pini{k, Zi„i) the 
initial conditions for N-body simulations are generated with 
a modified version of the GRAFIC2 package (Bertschinger 
2001). In this modified version the transfer function at the 
starting redshift is given as an external input. The start- 
ing redshifts z; are set to the time when the standard de- 
viation of the smallest density fluctuations resolved within 
the simulation box reaches 0.15 (the smallest scale resolved 
within the initial conditions is defined as twice the intra- 
particle dist ance). We used th e best-fit cosmological parame- 
ters fro miBovarskv et al.l (l2009a ). comparable with WMAP5 
results l|Komatsu et al.l 120091 '). Namelv: flm = 0.253, ~ 
0.747, n = 0.973, h = 0.72, and ag = 0.8. 

In our small box simulations (L = 20 h~^ Mpc ) we also 
include the effect of a non zero primordial velocity dispersion 
for WDM particles. The streaming velocities were generated 
using the corresponding expression for the average primor- 
dial vel ocity {v{z)) given e.g. in Section 4 in lBovarskv et al] 
(|2009al ): 
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Box size 
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part, mass 


force soft. 




[h-lMpc] 




[/i-1Mq] 


[h-^kpc] 


WDM-20 


20 


350^ 


1.31e7 


1.42 


WDM-45 


45 


350^ 


1.49e8 


3.21 


ml-20vel 


20 


350^ 


1.31c7 


1.42 


ml-20 


20 


350^ 


1.31e7 


1.42 


ml-45 


45 


350=^ 


1.49e8 


3.21 


ml-90 


90 


400=* 


7.99e8 


5.62 


ml-180 


180 


400=^ 


6.39e9 


11.25 


m5-20 


20 


350^ 


1.31e7 


1.42 


m5-45 


45 


350=* 


1.49e8 


3.21 


m5-90 


90 


400=* 


7.99e8 


5.62 


m5-180 


180 


400^ 


6.39e9 


11.25 



Table 1. N-body simulation parameters 
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{v(z)) = F„di 



3.15ir^(z) 

1 
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where z is the redshift of the beginning of simulation, m„dm 
is the mass of the WDM particle. We had kept the ampli- 
tude of primordial velocity constant and equal to ((3| and 
randomly chosen the direction of velocities of the particles 
to add them to the Zel'dovich velocities, generated with the 
GRAFIC2 package. 



3 NUMERICAL SIMULATIONS 

All simulations have been performed with pkdgrav, a tree 
code written by Joachim Stadel and Thomas Quinn (Stadel 
2001). The code uses spline kernel softening, for which the 
forces become completely Newtonian at 2 softening lengths. 
Individual time steps for each particle are chosen propor- 
tional to the square root of the softening length, e, over the 
acceleration, a: Ati — rjy^ e/oi. Throughout, we set r] — 0.2, 
and we keep the value of the softening length constant in 
comoving coordinates during each run. The physical values 
of e at z; = are listed in Table[T] Forces are computed using 
terms up to hexadecapole order and a node-opening angle 
6 which we change from 0.55 initially to 0.7 at z = 2. This 
allows a higher force accuracy when the mass distribution is 
nearly smooth and the relative force errors can be large. 

Table[T]lists all of the simulations used in this work. We 
have run simulations for several different box sizes, which 
allows us to probe halo masses covering the entire range 
10^°h-^ M© ^ M s; 10"h"^ Mq. For the smaU box simula- 
tions (L = 20 Mpc/h) we have, for each model, two different 
runs, with and without thermal velocities. 

The difference on small scales is lower in the evolved 
power spectrum than in the initial one (see fig. [T]). This is 
due to the non linearity of structure formation that is able 
to increase matter power on small scales, as shown by the 
regrowth of the power spectrum for k > 200 h/Mpc. As time 
goes by the scale at which non linearity dominates moves 
toward smaller values of fc, as shown by the comparison of 
the z = 2 and z — spectra. 

Fig. [2] shows the halo mass function for the m5 model 
at z = 0. Only haloes with viral masses larger than 
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Figure 2. Halo mass function for all three models at 2 = 
follow predictions of ACDM wi th the same cosmolo gical param- 
eters (green dashed line is the IWarren et alj ||2006| ) predictions 
for ACDM). Shown are the points for m5 model, the other two 
models (ml and WDM) give very similar results and are not 
shown. 



2 X 10^" Mq) are shown, and different symbols (and col- 
ors) refe£_to_dj5eren^box size simulations. The green line is 
the I Warren et all (j2006l ) prediction for a pure ACDM model 
with the same cosmological parameters as our CWDM mod- 
els. 

The I Warren et al] (|2006l ) prediction is a good representation 
of our data at high and low masses. The ml and WDM 
models present similar behaviours as the m5 one. 



3.1 Halo parameters 

In all of the simulations, dark matter haloes are identified 
using a spherical overdensity (SO) algorithm. We use a time 
varying virial densit y contrast det e rmine d using the fitting 
formula presented in iMainini et all (|2003l ). We include in the 
hal o catalogue all the h aloes with more than 500 particles 
(see lMaccio et all (j2008') for further details on our halo find- 
ing algorithm). For each SO halo in our sample we determine 
a set of parameters, including the virial mass and radius, 
the concentration parameter, the angular momentum, the 
spin parameter and axis ratios (shape). Below we briefly de- 
scribe how these parameters are define d and determined. A 
more detailed discussion can be found in Maccio et al.l l|2007l . 
|2008| ). Finally following iMaccio et al] |20oJ7'we~split our 
halo sample into unrelaxed and relaxed haloes. In the rest 
of the paper we will only discuss the properties of relaxed 
haloes. 
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Figure 3. Density profiles for CWDM and WDM models. The 
plot shows the change in the inner slope of the halo density profile 
for one of the haloes in the simulation with the box size L = 
20 Mph/h (as compared to the CDM - black line). Points and 
red line (the NEW fit) are for the m5 model. For ml and WDM 
models we show only the NFW fit (to make the plot less crowded). 
Black line sh ows the predicted profi l e for A CDM using the fitting 
formula from iMufioz- Cuartas etlH ||2010|'I . 



sity {pi) in 50 spherical shells, spaced equally in logarithmic 
radius. Errors on the density are computed from the Pois- 
son noise due to the finite number of particles in each mass 
shell. The resulting d ensity profile is fit with a NFW profile 
l|Navarro et al.ll 19971 ) Q: 



pir) 

Pc 



Sc 



{r/r,){l + r/n)2 



(5) 



During the fitting procedure we treat both rg and 5c as free 
parameters. Their values, and associated uncertainties, are 
obtained via a minimization procedure using the Leven- 
berg & Marquardt method. We define the r.m.s. of the fit 
as; 



Prms 



1 



In p„ 



(6) 



where pm is the fitted NFW density distribution. Finally, we 
define the concentration of the halo, Cvir = r^ir/ra, using the 
virial radius obtained from the SO algorithm, and we define 
the error on logc as (ar^/rs)/ In(lO), where is the fitting 
uncertainty on r^. 



3.1.1 Concentration parameter 

To compute the concentration of a halo we first determine 
its density profile. The halo centre is defined as the location 
of the most bound halo particle, and we compute the den- 



^ In this paper we do not resolve (due to resolution) and do not 
di scuss the inne r slop e of DM density profiles, as for example 
in iMaccio et alj ||2012| V Therefore, possible deviations from the 
NFW profile due to WDM effects, even appearance of a core, are 
not considered. 
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Figure 4. Mass and redshift dependence of the concentration parameter. Points with error-bars are CWDM si mulation (m5 and ml 
model s at left and right panel correspondingly) results, the straight lines are the fits from ACDM simulations of iMuiioz-Cuartas et al] 
Prom top to bottom: z = 0, z = 0.5, 2 = 1. One sees that ml model is "colder" (closer to the CDM) than m5 at all redshifts. 
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Figure 5. Same as fig. |4]for the WDM simulation 



3.1.2 Shape parameter 

Determining the shape of a three-dime nsional distributio n 
of particl es is a non-trivial ta sk (e.g., Ijing fc Sutdl20o3 ). 
Following lAUgood et al](|2006l ). we determine the shapes of 
our haloes starting from the inertia tensor. As a first step, we 
compute the halo's 3x3 inertia tensor using all the particles 
within the virial radius. Next, we diagonalize the inertia 
tensor and rotate the particle distribution according to the 
eigenvectors. In this new frame (in which the moment of 



inertia tensor is diagonal) the ratios s = a^/ai and p 
a-ijax (where ai ^ 02 ^ 03) are given by; 



03 

Hi 



03 

a2 



E '"if? 



(7) 



Next we again compute the inertia tensor, but this time 
only using the particles inside the ellipsoid defined by ai, a2, 
and a-i. When deforming the ellipsoidal volume of the halo, 
we keep the longest axis (ai) equal to the original radius of 
the spherical volume (rvir). We iterate this procedure until 
we converge to a stable set of axis ratios. 



4 THE CONCENTRATION MASS RELATION 

Figure |4] (left) shows the mass and redshift dependence of 
the concentration parameter for the m5 model. Results for 
2; = 0, 2 = 0.5 and 2 = 1 are shown from top to bottom. The 
error-bars on the points show the scatter in the concentra- 
tion mass relation, while the three li nes are the ACDM re- 
sults at different redshifts taken from 

Im ufioz-Cuartas et al.l 

(|2010l) 

At high masses the CWDM models basically agree 
with the pure cold dark matter predictions, even if there 
is an indication of a slighter higher concentration for the 
CWDM model for masses M ^ 10"/i"^ Mq, but still within 
the error-bars. The situation is different for massed below 
10^^/i~^Mq, where the CWDM model predicts lower con- 
centrations with respect to the ACDM. For our lowest mass 
bin (M « 10^°h'^ Mq) the concentration is lower by 45% 
and there is a clear indication of flattening of the con- 
centration mass relation, with an almost fiat relation for 
log(A/) < 10.5. This fiat tail of the cM relation at low masses 
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is already in place ai z = 1 and the difference between CDM 
and CWDM appears to be redshift independent. 

This results is not surprising. In CWDM models the 
formation of small haloes is delayed with respect to CDM 
due to the lack of power on small scales. The concentration 
parameter is related to the density of the univ erse at the 
time of the halo formation l|Wechsler et al.|[2002l ): since the 
density of the universe decreases with time, a later formation 
time implies a lower value for the concentration. 

The ml model behaves similarly to m5 as shown in 
right panel of figure |4] a slightly higher concentration at 
high masses and a less flatter cM relation at low masses. 
In this case the difference with ACDM is less pronounced 
and at the lowest mass scales (M ~ 10^^ Mq) it is less 
than 10%. The difference between ml and m5 is a direct 
consequence of the different power spectrum at small scales. 
As shown in Figure[T]the ml model has more power at scales 
fc Ri 1 — 5 h/Mpc, which correspond to mass scales of the 
order 1 - 5 x 10"/i"^ Mq. The higher power makes it more 
similar to LCDM and results in earlier formation time these 
halos and hence higher concentration. 

This is also confirmed by fig\5\ where the concentration- 
mass relation for the WDM simulation is presented. Results 
for the pure WDM model are very similar to the m5 model, 
confirming the relation between initial power spectrum (fig 
[TJ and halo concentration. 

4.1 Effects of thermal velocities 

For the small boxes (20 Mpc/h) we also run an additional 
simulation that included a thermal velocity component in 
the initial conditions. This thermal component is few percent 
of the initial velocity due to the potential field (according to 
the Zeldovich approximation), nevertheless it is important 
to test its effects (if any) on the halo internal structure. 

Figure [5] shows the one-to-one comparison of the me- 
dian Cvir at a given halo mass for the model ml, with and 
without primordial thermal velocities. As expected given the 
magnitude of the thermal velocity component, the effect is 
negligible. The model m5 presents the same behaviour. The 
thermal velocity component will play a role on smaller scales 
(for example for Milky Way satellites) and we plan to ad- 
dress this issue in a forthcoming paper. 

4.2 Effects on halo shape 

There is a well known relation between halo shape and 
mass with low mass haloes being more spherical than 
high mass ones (lAllgood et all l2006l : iMaccio et all l2008l : 
[m unoz- Cuartas et al.l boid '). This relation is usually ex- 
plained by assuming that the halo "sphericity" correlates 
with formation time, such that early halos have more time 
to virialize and hence reach a more equilibrated, less triaxial 
configuration. If this is the case we should expect to see a 
different trend between halo shape and mass in the CWDM 
simulation with respect to the ACDM. 

Figure [7| shows the redshift evolution of the halo shape, 
quantified via the minor to major axis ratio: s = a3/al, for 
the ml and m5 models (left and right panels in Fig.[7]corre- 
spondingly). On average it looks like that at 2: = low mass 
haloes are slightly more spherical in CWDM than in the 
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Figure 6. Effects of thermal velocities on the concentration pa- 
rameter in the 20 Mpc box. Different symbols refer to differ- 
ent redshifts. 

CDM, while the trend is somehow reversed at high masses. 
This is consistent with the results on the concentration and 
confirms that the halo formation time is a key ingredient in 
defining the structural halo properties. 



5 SUMMARY AND DISCUSSION 

Interest in warm dark matter models has increased in 
the last years, mainly motivated by the possible need 
for cored (or very low concentration) dark matter halo 
density p rofiles for dwarf galaxies i n the local group 
(see e.g IWalker fc Penarrubial I2OIII : lAmorisco fc Evana 
l201ll ). or b y solution of vario u s CDM "ove r-abundance 



Droblems" (iKlvDin et all Il999l: 


Moore et all IT999l~). see 


e.e. jMaccio & FontanotI 201C 


: Strigari et al. 


201c ; 


Bovlan-Kolchin et al.l I2OIII: iPaoasterEis et al.l 


2011; 


Lovell et al. 20111). 







While it is still under debate what is the effects of 
bary ons on the dark matter distribution on such sma ll scales 
(e.g. iPontzen fc Governatdl20l"ll : iMaccio et al ]|2012l and ref- 
erences therein); it is worth studying through numerical sim- 
ulations the effects of a warm component on the inner struc- 
ture of DM haloes. It had been repeatedly argued that "pure 
warm" dark matter models (with a cutoff-like suppression 
of the matter power spectrum) could not be a solution to 
the core-cusp problem as the required free-streaming would 
be in a stark contradiction with Lyman- a bounds (see e.g. 
IStrigari et al]|2006l ; Ide Narav et"aLll2009l ) . 

In this paper we have performed the first N-body sim- 
ulations for cold plus warm dark matter models (CWDM). 
In this class of models the dark matter particle candidate 
(e.g. sterile neutrinos, gravitino, axions) can be produced 
via two co-existing mechanisms and therefore the primor- 
dial spectrum is a superposition of a cold and a warm com- 
pone nt, plus a complicated non-thermal velocity spectrum 
(e.g. iBovarskv et~aLll2009bh 
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Figure 7. Mass and time evolution o f tlie shape of dark matter i ialoes quantified via tiie s = a3/al ratio. The black dashed lines are 
the results for the ACDM model from lMuiioz-Cuartas et al. Left panel - ml model, right panel - m5 model 



We have mainly focused our attention on the inner 
structure of dark matter haloes in CWDM models, and its 
evolution with redshift. We have adopted the commo nly 
used concentration parameter (e.g. iMaccio et al.l l2008h to 
parameterize the modification in CWDM models with re- 
spect to the standard Cold Dark Matter (CDM) model. 

We show that in our models the number of halos re- 
mains the same as in ACDM down to the smallest scales 
resolved in the simulations (~ 10^" Mq). At the same time, 
keeping the free-steaming scale the same, the variation of 
WDM fraction is able to reduce the concentration parameter 
on mass scales as high as M « 10^^h~^ M0, two/three order 
of magnitude above the free streaming mass of the model. 
On dwarf galaxies mass scale(M ~ 10^^ Mq) the concen- 
tration parameter is half of the value predicted by ACDM 
and we see a clear sign of a strong flattening of the concen- 
tration mass relation. As a consequence we expect an even 
stronger concentration reduction on lower mass scales (be- 
low our resolution limits), mass scales directly explored by 
the Milky Way satellites. This decrease of the concentration 
of such large halos may explain why the "cuspy" matter dis- 
tributions are not supported by ob servations of the rotation 
curves of spira. 1 galaxies (see e.g. Salucci fc Burkert I2OO0I : 
lOh et al.l I2OO8I : ISpano et al.l boosl : Ide Narav et all |2009^ ■ 
thus resolving one of the major challenges for the CDM cos- 
mological model. 



In lMaccio et al.l l|2012h it was shown that the primordial 
velocities required to produce the cores of observable sizes 
in pure WDM cosmologies are so large that too few galaxies 
would be formed. This is because in pure warm dark matter 
models the phase-space density, defining the size of the core 
and suppression in the halo number density are defined by 
the same parameter - the average velocity of particles. Our 
results show that in the CWDM models the presence of the 
second parameter, F„dm, allows to control the "warmness" 
of dark matter, separating the scales of modification of the 



matter power spectrum (halo number density) and that of 
modification of the density profile (size of the core). 

Our analysis also suggests a slightly higher concentra- 
tion for massive haloes with respect to CDM (but still within 
the error). This higher central density would help in recon- 
ciling DM based predictions with, results fr om X-Ray and 
gravitationa l lensing (see e.g. [ Comerford fc Nataraian. 20071 
but see also lKlvpin et aLll2011 ). 

We stress that unlike the recent work of Im accio et al.l 
(I2OI2I ) this effect is not related to the finite phase-space den- 
sity (the phase-space density of particles in our simulations 
is the same as in pure CDM ones). We ascribe the lower con- 
centration at small masses to a shift in the halo formation 
time, that results in a delayed formation for CWDM models 
with respec t to CDM, as already f ound in previous studies 
(|Eke et al.l^OOll : lLovell et allboill ). This delayed formation 
time has also an effect on the halo shape, with CWDM ha- 
los being, on average, more triajcial than CDM ones. We did 
not find any effect on the halo spin parameter distribution, 
which looks almost identical in all models. 

We also tested the effect of explicitly including thermal 
velocities in the N-body initial conditions. We did not see 
any visible change in the concentration mass relation on the 
mass scales we are able to resolve. We plan to study in more 
details the importance of this thermal velocity component 
in a forthcoming paper, by employing "zoomed" high reso- 
lution simulations of single dark matter haloes. 

Along with the Lyman-a forest method, the weak lens- 
ing surveys can be used to probe further clustering prop- 
erties of dark matter particles as sub-galactic scales, as 
the next generation of these surveys (such as e.g. KiDS, 
LSST, WFIRST, Euclid) will be able to measure the matter 
power spectrum at scales down to 1 — 1 h/Mpc with a few 
percen t accuracy. iMarkovic et al.l (|201l[ ) ; ISmith fc Markovi3 
j201ll ) argued that the next generation of lensing surveys can 
provide sensitivity, compatible with the existing Lyman-g 
bounds (|Viel et al.ll2"006l : ISeliak et al.ll2006l : iBovarskv et all 
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l2009al )). As in the case of the Lyman-a forest method the 
main challenge for the weak lensing is to properly take into 
account baryonic effects on matter power spectrum. The 
suppression of power spectrum due to primordial dark mat- 
ter velocities can be extremely challenging to disentangle 
from the modificat ion of the matter power spectrum due to 
baryo nic feedback (jSemboloni et al1l201ll : Ivan Daalen et al.l 
I2OIII ). Finally, the modified concentration mass relation 
can be probed with the weak lensing sur veys (see e.g. 
iMandelbaum et al.ll2008l : iKing fc Meadlboill ) if their sensi- 
tivity can be pushed to halo masses below roughly 

While the observational difference between pure cold 
and mixed cosmologies is not drastic, clarifying this issue 
would have profound impact on the fundamental physics 
questions. Cold-|-Warm dark matter models are still in an 
infant state and an effort comparable with that, invested in 
theoretical investigation of structure formation of pure CDM 
may be needed before this question will be finally settled. 
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